Implementacja w MATLABie problemu transportu ciepła
Poniżej przedstawiam kod MATLABa obliczający problem transportu ciepła na obszarze w kształcie litery L.
Pobierz kod
lub zob. Załącznik 3.
Kod może zostać uruchomiony w darmowym środowisku Octave.
Kod uuchamia się otwierając go w Octave oraz wpisując komendę
\( heat' \)
Po chwili obliczeń kod otwiera dodatkowe okienko i rysuje w nim rozwiązanie.
W linii 641 \( knot = [0, 0, 0, 1, 2, 2, 2]; \)
znajduje się zapisany wektor węzłów [0 0 0 1 2 2 2] który używany będzie na każdym z trzech grup elementów (patchów) stosowanych do rozwiązania problemu transportu ciepła.
Treść zadania:
Proszę zmodyfikować kod transportu ciepła, tak żeby używał on wielomianów Lagrange'a drugiego stopnia, i żeby każdy patch elementów miał 5 elementów
Autorzy kodów w MATLABie: Marcin Łoś i Maciej Woźniak.